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Abstract 

We investigate the pairing property of nuclear matter with Relativistic 

Hartree-Bogoliubov(RHB) approach. Recently, the RHB approach has been 

widely applied to nuclear matter and finite nuclei. We have extended the 

RHB approach to be able to include non-linear coupling terms of mesons. In 

this paper we apply it to nuclear matter and observe the effect of non-linear 

terms on pairing gaps. 
21.60.-n, 21.65. +f 



Typeset using REVTeX 



*e-mail: satoru@postman.riken.go.jp 
^e-mail: sumi@la.numazu-ct.ac.jp 
^e-mail: toki@rcnp.osaka-u.ac.jp 



1 



I. INTRODUCTION 



The development of secondary beam to produce unstable nuclei changed the common 
sense of nuclear physics, which had been accumulated over the years from the study of 
stable nuclei. The discovery of neutron halo and skin is one of the most famous ones. These 
observations changed the formula of the nuclear radius being proportional to A 1 / 3 . New 
facilities are now planned and being constructed to get more unstable nuclei including those 
near the drip line. These new facilities would definitely bring us with more information and 
new phenomena for those unstable nuclei. 

In the theory side, various attempts were made to understand diverse phenomena with 
the unified frameworks. One of the powerful framework is the relativistic mean field the- 
ory (RMF) It is based on the relativistic many-body theory and is applied to many 
observables with great success. 

In the RMF theory, it is very important to include the pairing interaction on the same 
footing as the particle-hole correlations. Usually, it is implemented by the so-called BCS- 
type theory as in other mean field theory. In the BCS-type theory the pairing interaction 
is implemented by hand outside of the relativistic theory. However, we would like that the 
pairing interaction is described in the relativistic many-body theoretical manner within the 
RMF theory. One possibility toward this attempt is the relativistic Hartree-Fock-Bogoliubov 
theory(RHFB), which was formulated by Kucharek et al f|. Recently, the Hartree-Fock- 
Bogoliubov(HFB) theory is widely used in both the non-relativistic and relativistic mean 
field calculations. It can treat the usual mean field and the pairing field simultaneously. 
The HFB theory is applied to nuclear matter and finite nuclei, relativistically PHT^I and 
non-relativistically [|T~3| |3~5]| with many interesting results. 

However, until now the effect of the non-linear self-coupling terms of mesons are not 
considered in the pairing channel in the RHB theory. It is known that those non-linear 
terms are needed to reproduce the properties of nuclear matter and finite nuclei in the RMF 
theory |]I6|-p^] . We can imagine that those terms have some effect on the pairing correlations. 
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Therefore, in this paper we would like to extend the relativistic Hartree-Fock-Bogoliubov 
theory to include the effect of the non-linear self coupling terms of mesons and apply it to 
nuclear matter to see the effect on the pairing gap A. 

In Section [TI], we develop the formalism of the RHFB theory using the non-linear rela- 
tivistic Lagrangian under the harmonic approximation to handle the non-linear terms. In 
Section |T|, this equation is solved with the parameter sets of TM1 [18| and NL1 | 19fl . Since 
these parameter sets were obtained in the RMF formulation, the Fock term contribution is 
dropped in this paper. We study the form factor dependence on the pairing properties of 
the nuclear matter. Section |IV| is devoted to the summary of this paper. 



II. FORMALISM 



In the RMF theory we use the Lagrangian density p|-|^JT8[], 



C(x) — if> (x) 



i^d? — M — g a a — g^l^ - # p 7 m t> 



fj,a 



1-T 3 



if) (x) 



+ X -d^ (x) d»a (x) -U{a)- \v, u V^ + ^m 2 w u^ + ±c 3 (lu^) 



(1) 



In this Lagrangian if) is a Dirac spinor of a nucleon. We introduce three types of meson fields 
a, u) and p. a is the isoscalar-scalar field, which mediates medium range attraction between 
two nucleons. u is the isoscalar-vector field, which produces short range repulsion, p is the 
isovector-vector field, which gives isovector correlation among nucleons in nuclei. A is the 
photon field, which produces the electromagnetic field. Here, U (a) is self coupling terms of 
the cr meson, 



U (<r) = \™W + \92^ + ^g 3 cr 4 . 



(2) 



Boguta et al. investigated it extensively [IBj and found that it is needed to reproduce the 
nuclear matter and finite nucleus properties. We also include a non-linear coupling term of 



co meson |17| , P^8| , which was investigated by Sugahara et al., motivated by the recent success 



of the relativistic Briickner-Hartree-Fock(RBHF) theory ||20|| . They found that it is needed 



to reproduce the RBHF result for the nuclear matter quantitatively up to the high density 
region. From this Lagrangian, we can derive a Dirac equation for nucleons and Klein-Gordon 
equations for bosons with the Euler-Lagrange equation. In the RMF theory we replace the 
meson fields with their classical mean values, 

a — > (a) etc. (3) 

However, in this approximation the terms which lead to the pairing interaction do not 
appear. These terms come from the terms like ij)* ip' 'ipip . To obtain these terms we need to 
take the higher order terms and fluctuation terms of the meson fields ||. Here we explain 
our method which is somewhat different from Kucharek et al. For brevity, we write down 
only the a meson terms. The extension to include other mesons is straightforward. We write 
the quantized field as the sum of its classical field and its fluctuation, 



a = a + a 1 '. 



(4) 



Here, the classical mean value a is obtained from the mean field equation, 

g a + Da + U\a) = 0. (5) 
By substituting Eq. (|4]) into the relativistic Lagrangian ([!]), we get, 
C (x) (x) {h^ -M*}i(j (x) 



- g a i> (x) i> (x) a' (x) + ^d^d^x) - U {a) 



+ d^a' (x) d^a (x) - U' (a) a' (x) 

+ U,a' (x) d»a' (x) - l -U" (a) a' 2 (x) + 0(a' 3 ). (6) 

0(cr' 3 ) indicates the terms of more than a' 3 orders. We neglect the higher order terms, 
0(a' 3 ), i.e., the harmonic approximation. In the above equation M* is the effective mass of 
a nucleon, 

M* = M + g a a. (7) 



In Eq. (||), we introduce the new mass parameter \i for the a' field, 

= U"{a) = m 2 a + 2g 2 a + 3g 3 a 2 . (8) 

Using this, the propagator for a' meson field can be written, 

/dk 4 e ~ ik ( x ~y) 

It is noted that the mass term in the a' meson propagator is /i, not m a . Using this, we can 
write down a' field, 

a'(x) = J d 4 yD(x - y) {g a $(y)ij>(y) + U<t(y) + U'{a)) . (10) 

Substituting this expression into Lagrangian (||), we get, 

x D(a: - y) {gMy)1>(y) + □<t(</) + £/'(*)} . (11) 

We may consider this effective Lagrangian as obtained by integrating out the fluctuation 
field a'. For the interaction part we use the Gor'kov method fSffil, which is used in the 
theory of superconductor , 



d 4 xCi n t = — J j d 4 xd 4 yD(x — y) {g a if)(x)if)(x) + Da(x) + U'(cr(x))} 

x {g„ (HyMy)) + □%) + U'(a(y))} 

d 4 xd 4 yD(x - y)i>{x) (ip(x)i)(y)) ip(y) 



+ 2& II d 4 xd?yD{.c - </)<•(.<■)<'(</) ■(< m ./■)(-(//)) 



+ -gl J J d 4 xd 4 yD(x - y) {${x)${y)) ^{x)^{y). (12) 

In the above equation, the first term vanishes because of the mean field equation (|^). The 
second term corresponds to the Fock term and the last two terms are responsible for the 
pairing interaction. Finally, we get a modified effective Lagrangian of the form, 

^ e// = £-MF + Cint- (13) 



In this Lagrangian Cmf corresponds to the Lagrangian of the RMF theory, 



MF 



$(x) (iy„d» — M — g a a{x)) ^) + \d.a{x)d^a{x) - U(a), 



(14) 



with which we can reproduce the RMF theory. £i nt is given in eq. (|T2|). 

Using the Euler-Lagrange equation to this Lagrangian, with respect to if) and ip we get 
the Dirac equations for the nucleon field ij) and ip, 



(i 7 ^ - M*) ab Mx) =gl J d%D( x - y) (M*)Mv))Mv) 

-gl d 4 yD(x-y) (if) a (x)il> b (y)) $ b (y), (15a) 



- M*) ba Mx) =g 2 a J d A yD{x - y) (My)Mx)) My) 

-gj [ d A yD{x - y) (Mx)My)) My)- (15b) 



With these equations we get the relativistic Hartree-Fock-Bogoliubov equation as written in 
Kucharek et al. 



\ 



h-X A 
-A* -h* + Xj 
Here, h is the Dirac Hamiltonian and A is pairing field, 



(16) 



hob (x, y) 



6{x-y) 



A ab (x,y) = J2 f d 3 z ( 7 ^ (x,z)) 



abed 



Kcd (y,z) 



(17a) 
(17b) 



Here, we have included the chemical potential A which is fixed by < A^t^A >= N with 
N being the particle number. In the above equations, p and k are the normal density and 
the pairing density. 



Pab (x,y) = (A\if>l (y)ip a (x)\A), 
Kab (x,y) = {A\ip b (y) ipa (x)\A). 



(18a) 
(18b) 



V^'s are the meson potentials, which are introduced in the same procedure as the case of the 
a meson: 

VU*. V) = Si (iT*). (7°r*) M / (g, (* = •>, <-, P, A) . (19) 

is the nucleon-nucleon-meson vertex, for example, T u = 7^. The only difference between 
Kucharek et al. and ours is that in our formalism the masses of mesons in p-p channels are 
modified by non-linear terms of mesons, (see Eq (§).) 



III. RESULT 

In this section we apply the formalism derived in the previous section to nuclear matter. 
Before we mention the result, we explain the approximations that we use in this work. First, 
we neglect the Dirac sea as usually done in the RMF theory. Second, we neglect the Fock 
term appeared in Eq. (|12|) as usual because our aim is to check the effect of the non-linear 
terms of mesons on the pairing channel. Third, we take into account only the 1 <So pairing, 
which is the most important one in both nuclear matter and finite nuclei. We do not treat 
proton-neutron pairing. With these approximations, we deduce from Eq. fll6|) the BCS-type 
gap equation ||, 

1 r n. \ a w 



A ^ = -^ %,(M~j= Pdp. (20) 

Here, A(/c) is the pairing gap in the nuclear matter at momentum k. A is the chemical 
potential. v pp (k,p) is the interaction term of the pairing channel. As for the precise form of 
v P p(k,p), see Kucharek et al H. Furthermore, we introduce the form factor for each vertex, 

A 2 _ 777 2 

F a (q 2 )= °_ f (a = a,^andp) (21) 

a 9 

to cut out the effect of the high momentum region. Without the cut-off, we get unreasonably 
large pairing gap |||8|]. In this work we use the cut-off parameter A, 1800MeV for a, 1500MeV 
for u) and 1300MeV for p. These form factors correspond to those of Bonn C [2(J. We note 



that there is large cut-off dependence of the value of the pairing gap. However, we find that 
the overall tendency is not changed. The actual form of v pp is plotted in Fig [TJ. In this figure, 
we also plot the components of the a and cu meson exchanges. We use the TM1 parameter 
set . From this figure, we see that the large attraction of a and the large repulsion of uj 



cancel out each other and the relatively week interaction in the lower momentum region is 
realized. It is one of the characteristic feature of the RMF theory. 

In Fig. H] we show the pairing gap at the Fermi momentum kp fm _1 as a function of 
the density measured in the Fermi momentum. We use the parameter set TM1, which has 
both a non-linear coupling and uo non-linear coupling terms ||18|| . In this figure, "NL" (solid 
line) denotes the case when we use the formalism in the previous section to include the 
non-linear term effect on the pairing channel and "L" (dashed line) means the case with the 
use of the raw mass m a , m u for the mass parameters in the propagators of mesons instead 
of \i a and (see Eq. (0)). From this figure we see that the pairing gap is reduced by a few 
MeV by including the non-linear term of mesons. This effect is especially large around kp 
= 1.0 fm -1 . Only the difference between two cases is the mass terms in the propagator of 
mesons. In our formalism, by including the non-linear term effect the masses of mesons in 
the propagator are modified with the medium effect. To see it, we plot (solid line) and fi u 
(dashed line) in Fig. || At lower density both \x a and /i^ are not much different from their 
original values, m a and m^, because the densities are low, medium effect is not large in this 
region. However, at higher density, they are larger than their original values. These results 
come from the non-linear terms in the expressions of and ^ (Eq. (§)). Those effects 
are more significant at higher density. Because the mass terms in the propagator are more 
massive at higher density region, the interaction in the pairing channel is reduced. To see 
that in Fig. |1| we compare the value of v pp at the density kp = 0.78 fm -1 , which corresponds 
to the maximum point of A(kp) without the non-linear term effect in Fig. |2|. We see that by 
introducing the non-linear term effect, v pp is reduced. Hence, pairing gaps are reduced. Of 
course, those effects depend on the value of coupling constant. However, at normal density 
the general forms of non-linear terms are not so different between the parameter sets which 
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are widely used. So, we can find the same tendency; namely, the reduction of pairing gap 
with each parameter set at moderate density. 

We study the effect of the form factor by changing the cut-off mass slightly. With the 
use of A CT =2200MeV, A w =1900MeV, A p =1700MeV, the pairing gap increases largely as 
shown in Fig The peak position, however, does not change much from the previous choice 
of the parameters. These cut-off masses in the corresponding form factors are related with 
the structure of hadrons in the nuclear medium. It is certainly an interesting future subject 
to fix these values on theoretical ground. 

For comparison we show A(kp) for the case of the NL1 |19] parameter set in Fig. |5]. NL1 
has the non-linear term of only a meson. It does not have those of u meson. We can see 
that the same effect appears as in the case with TM1. However, at higher density A(kp) 
get abnormally large values. (It is not shown in the figure.) This is due to the minus sign of 
the coupling constant of the quadratic self-coupling term. Because of this negative sign, at 
higher density fi a is getting smaller and smaller with density. Hence, the pairing interaction 
becomes very large and the pairing gap becomes very large at high densities. In the TM1 
case, because g 3 has positive sign this effect does not occur. Nevertheless, we can find the 
same tendency in the moderate mass region. 



IV. CONCLUSION 

We have developed a new framework to incorporate non-linear terms of mesons in pairing 
channel in the RHFB theory. We have applied the RHB theory to nuclear matter to study 
the effect of non-linear terms on the p-p channel. We have found that the non-linear terms 
reduce A at Fermi momentum especially for normal density region. It is caused by non- 
linear term effect on the mass parameter in the meson propagator. The pairing interaction 
is reduced due to the non-linear coupling terms. This effect is neglected in other works. We 
have found that it has a non-negligible effect. It is interesting to apply it to asymmetric 
nuclear matter and finite nuclei. 
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In this paper we have used the parameter sets determined in the RMF framework; i.e., 
TM1 and NL1. As seen in the formalism, the pairing correlations in nuclear matter are 
caused by the fluctuations of the meson fields. Hence, it is very natural to take the Fock 
term when we consider the pairing correlations. It then needs an elaborate work to complete 
our program as to first fix the parameters of the relativistic Lagrangian in the Hartree-Fock 
framework from various magic number nuclei and then to calculate the pairing properties 
of nuclear matter. We start this program to be worked out in due time. 
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FIGURES 

FIG. 1. The interaction kernel in the gap equation v pp (p, k) at p = 0.78 fm _1 for the case of 
the Fermi momentum, kp = 0.78 fm _1 . The horizontal axis is k in fin -1 . "NL" (solid line) denotes 
the case of non-linear terms of mesons being incorporated and "L" (dashed line) denote the case 
without them. "Total" corresponds to the sum of the contributions from a, u, and p mesons, where 
the curves with a and ui denote the contributions from a and u mesons, respectively. 

FIG. 2. Pairing gap A at the Fermi momentum &p as a function of the Fermi momentum. 
The TM1 parameter set is used. "NL" (solid line) corresponds to the case in which the non-linear 
term of mesons are implemented. "L" (dashed line) corresponds to the case without them. 

FIG. 3. The mass parameters of a and u) mesons in the propagators, \i a (solid line) and n u 
(dashed line). The horizontal axis is the Fermi momentum &f fm _1 . The increases of masses due 
to the non-linear terms of mesons are observed at higher Fermi momentum. 

FIG. 4. Pairing gap A at the Fermi momentum /cp as a function of Fermi momentum. The 
TM1 parameter set is used. Dashed line corresponds to the cut-off parameters, A CT =2200MeV, 
A w =1900MeV, A p =1700MeV. Solid one corresponds to the cut-off parameters, A CT =1800MeV, 
A^=1500MeV, A p =1300MeV. 

FIG. 5. Pairing gaps A at the Fermi momentum fcp as a function of Fermi momentum for 
the case of the NL1 parameter set. The solid curve corresponds to the case with non-linear term, 
whereas the dashed curve without the non-linear term. 
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